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Abstract 



In the simplest model of magnetic pairing, the transition temperature to 
the superconducting state depends on the dynamical susceptibility x(q, w). 
We discuss how T c is affected by different momentum and frequency parts 
of x(q ; ^) for nearly antiferromagnetic and nearly ferromagnetic metals in 
two dimensions. While in the case of phonon-mediated superconductivity 
any addition of spectral weight to a 2 F(uj) at u> > leads to an increase in 
T c , we find that adding magnetic spectral weight at any momentum q and 
low frequencies ([0 : 3T C ] and [0 : (5 — 9)T C ] for nearly antiferromagnetic 
and ferromagnetic metals respectively) leads to a suppression of T c . The 
most effective frequency and momentum range consists of large momenta 
q ~ (-7T, 7r) and frequencies around 10T C for nearly antiferromagnetic metals 
and small momenta q ~ and frequencies of approximately (13 — 22)T C for 
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nearly ferromagnetic metals. 
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I. INTRODUCTION 



In principle there are over one million ternary and over one hundred million quaternary 
crystalline materials. Even the binary compounds number well over eleven thousand. Only 
a tiny percentage of these materials have been synthesised. Clearly an exhaustive search 
of the periodic table for superconductivity is out of the question. However, there exist in 
nature a wide range of compounds which are very close to a magnetic instability, and in 
quite a number of cases, an anisotropic superconducting phase is experimentally found on 
the border of magnetic order. But this type of superconductivity is sometimes confined to 
very small regions of the phase diagram, making it difficult to detect. The heavy- fermion 
superconductor UG62 1 is a case in point. An intuitive understanding of the trends in T c not 
only provides a test of our models of anisotropic superconductivity but can also guide the 
experimental search for new anisotropic superconductors. 

One of the most extensively investigated models of anisotropic superconductivity is based 
on a magnetic interaction, in which quasiparticles interact via the exchange of ferromagnetic 
or antiferromagnetic spin-fluctuations. An understanding of the properties of this model is 
essential in assessing its ability to describe metallic systems on the border of magnetic 
long-range order. The role of dimensionality, whether one is close to a ferromagnetic or 
antiferromagnetic instability and the sensitivity of T c to model parameters have been inves- 
tigated within a mean-field framework. 2 ' 4 ' 3 In these calculations, the transition temperature 
reflects an integral property of the magnetic excitation spectrum and does not tell us how 
much different parts of %(q, uj) contribute to the answer. 

For the traditional isotropic superconductors, the transition temperature T c depends 
upon the spectral function a 2 F{uj) which characterises the phonon-mediated electron- 
electron interaction. The sensitivity of the transition temperature to variations in the spec- 
tral distribution a 2 F{uj) was first investigated by Bergmann and Rainer 5 . They found that 
all frequency regions of a 2 F(uj) yield a positive contribution to T c . Frequencies uj 3> T c 
and w < T c contribute little (nothing in the limit uj — > 0) while frequencies around 2ttT c 
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contribute the most. 

Later Millis, Sachdev and Varma 8 extended the Bergmann-Rainer analysis of the 
Eliashberg equations to the case of magnetically mediated singlet superconductivity They 
assumed that the dynamical magnetic susceptibility could be expressed as a momentum- 
dependent factor times a frequency- dependent factor and that the spin fluctuations could 
be described by an Einstein model. In their model the sensitivity of T c to changes in the 
magnetic spectral weight did not depend on the paramagnon momentum q and made the 
computational analysis of the problem very similar to that of the phonon case. They showed 
that, in contrast to the latter problem, there is a crossover frequency u cross such that adding 
magnetic spectral weight at frequencies uj < uJ cross led to a reduction in T c . 

More recently, Monthoux and Scalapino 6 carried out this analysis for the fluctuation- 
exchange approximation to the two-dimensional Hubbard model. For this more realistic mag- 
netic pairing interaction, which is very non-local in space, they generalized the Bergmann- 
Rainer approach in order to study how sensitive the transition temperature to the d x 2_ y 2 
superconducting state is to infinitesimal changes in the spectral weight at frequency uj and 
momentum q. In agreement with the results of Millis et al 8 , they found that adding spectral 
weight at very low frequencies and for any momentum q led to a reduction in T c . They also 
found that there was a region in momentum space near q = where any small addition 
of spectral weight at any frequency also led to a reduction in T c . They mapped out the 
region in q and uj space that gave a positive contribution to T c , essentially wavevectors near 
q = (tt, 7r) and frequencies uj > T c . 

Here we report the results of a similar analysis using a parametrization of the effective 
interaction arising from the exchange of magnetic fluctuations. We consider systems on 
the border of antiferromagnetism and ferromagnetism. We study changes in T c brought 
about by small variations in the spectral weight at some wavevector q and frequency uj. 
We map out the regions of wavevector and frequency space where a small addition of spin- 
fluctuation spectral weight yields an enhancement or suppression of T c . This will give a 
more detailed understanding of the trends in T c that one calculates from the model. We 
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shall contrast these findings with those obtained for the conventional phonon-induced pairing 
mechanism 5 to gain insights into the similarities and differences between phonon-mediated 
and magnetically mediated superconductivity. 



II. MODEL 

We consider quasiparticles on a two-dimensional square lattice and postulate the following 
effective action 



S *ff = Y,J Q dT ^pa ( r ) (dr + Cp - flj % a (t) 

jf /, d/x(q,r-/)s(q,r).s(-qy), (1) 



where ip pa . and ip pa are Grassman variables and N is the total number of allowed wavevectors 
in the Brillouin Zone. The spin density is given by 

s(q, r) = ]T ^ P +qa( r )<W'P7( r )> ( 2 ) 

pcry 

where cr is the vector whose components are the three Pauli spin matrices. 
The dispersion relation is 

e p = — 2t(cos(p x ) + cosfjOy)) — At' cos(p x ) cos(p y ), t' < 0.5t, (3) 

where for simplicity we use units in which the lattice spacing is unity. For comparison with 
earlier work on this model 2 ' 4 we set t' = 0.45t and adopt the value n = 1.1 for the band 
filling. 

We parametrize the retarded generalized magnetic susceptibility as: 

K 2 _|_ q 2 _ i[u/r)(q)\ 

where k and kq are the inverse correlation lengths with and without strong magnetic corre- 
lations respectively. The correlation length is related to the pressure applied to the sample 
with k 2 = coinciding with the quantum critical point. Let 



q± = 4 ± 2(cos(g x ) + cos(q y )). 



(5) 



For antiferromagnetic correlations the parameters q and r] are 



-2 -2 



v(q) = T sfQ-, 



where T s f is a characteristic spin-fluctuation temperature. In the ferromagnetic case 



-2 -2 

q =Q- 



(6) 
(7) 



(8) 
(9) 



We obtain the total spin-fluctuation propagator on the imaginary axis xi^^n) via the 
spectral representation 

r+oo ^ u Imx(q, cu) 



X(q, %v n ) = - — 

J-oo 7T 



2Z/„ — C<J 



(10) 



To get x(Q)^n) to decay as 1/i^ as z/„ — > oo, as it should, we introduce a cutoff oj cut and 
take Imx(q, u) = for u > u cut . A natural choice for the cutoff is u) cut = ii(q) k o- 

Using the effective action in Eq. (1) and the dynamical susceptibility in Eq. (10), the 
two-dimensional mean-field Eliashberg equations for the transition temperature T c in the 
Matsubara representation reduce to 



T 



E(p, iuj n ) = g 2 — ^2^2x{p- k, iu n - iQ n )G{k, iQ n ) 



G(p, iu r , 
A(T)$(p,^ r , 



iu n - (e p - n) - E(p,iw„) 



(11) 
(12) 



<7 2 /3 



-9 Z 



T 

]y 1212 X(P - k > iuJ n - i^n) 
^Ifi k 



x \G(k,iQ n )\ 2 $(k,iuj n ), 



(13) 



where 



A(T) = 1 =>- T = T c . 
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(14) 



S(p,ic<j n ), G(p,iu n ) and $(p, iuj n ) are the Fourier components of the quasiparticle self- 
energy, the one-particle Green's function and the anomalous self-energy respectively. In 
Eq. (13) the prefactor — g 2 is for singlet pairing while the prefactor g 2 /3 is for triplet pairing. 

We find an instability to a d-wave gap function with <E>(p, iu n ) transforming as cos{p x ) — 
cos(p y ) in the nearly antiferromagnetic case and an instability to a p-wave gap function with 
$(p, iuj n ) transforming as sin(p x ) or sin(p y ) in the nearly ferromagnetic case. 

The momentum convolutions in Eqs. (11) and (13) were evaluated with the aid of a fast- 
Fourier-transform algorithm on a 128 x 128 lattice. The corresponding frequency sums were 
carried out using the renormalisation group technique of Pao and Bickers 7 . Between 8 and 
16 Matsubara frequencies are kept at each stage of the renormalization group procedure. 
The renormalization procedure is started at a temperature T = OAt and the frequency sum 
cut-off used is Q c ~ 20t. The renormalization procedure restricts us to discrete tempera- 
tures so that the point at which the condition in Eq. (14) is met must be determined by 
interpolation. The discrete temperatures were sufficiently close that a linear interpolation 
was adequate. The renormalization procedure afforded us considerable savings in computer 
time and storage requirements. Because of this we were able to carry out a thorough analysis 
of the dependence of our results on the spin-fluctuation coupling parameter g 2 Xo/t an d the 
inverse correlation length parameter k 2 . 

To investigate how strongly the transition temperature is influenced by various frequency 
and momentum parts of the paramagnon spectral function, we add an infinitesimal amount 
of spectral weight at frequencies ujq > and —ujq and wavevector q* and numerically cal- 
culate the change in T c . More specifically, the paramagnon spectral weight is changed from 
Imx(q, u) to Im(x(q, u) + <5x(q, u)) with 

N 

lm5x{q,uj) = r]7i{5(u - u ) - 8{u + w ) } —— ^ 5^* , (15) 

q * q* 

where 77 is a positive infinitesimal dimensionless parameter. The sum over q* includes q* 
and those wavevectors in the Brillouin zone related to it by the symmetry operations of the 
lattice. iV q * is the number of such wavevectors including q* itself. iV is the total number of 
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allowed wavevectors in the Brillouin zone. 

In the Matsubara representation, the addition of this infinitesimal amount of spec- 
tral weight corresponds to changing the effective interaction V e ff(q, iv n ) = <? 2 x(q, i^n) i n 
Eqs. (If) and (13) to V eff + 5V eff where 

XT/ f ■ \ V 2 2uJ o N „ 

5V eff (q,iv n ) = —g Xo , 2^<W> ( 16 ) 

Xo ^0 + V n iV Q* q* 

The total effective interaction thus depends on the parameters k, k , g 2 Xo/t as we U as 
the ratio i]/xo, and therefore so does T c . For small 77, one has 

n dT 

T c (rj/ X 0: q*, w„, • • •) = T c (0, q*, cu , . . .) + — W^(0, q*, uj , ■ ■ ■) (17) 



where . . . denote all the other parameters in the problem. The quantity 

(18) 



AT c (q*, Wo ) = ^(q* )W o) 



r?=0 

is a measure of how sensitive T c is to an infinitesimal change in the value of Imx at (q*, cu ). 
We calculate this derivative using the finite-difference estimate 

dT, 



-* x _ T c (Vxo, q*, u ) - T c (0, q*, <Jq) 
<Wxo) »>= WXo 

The value of the parameter rj/xo must be chosen small enough, such that the function 

T c (r)/xo) is approximately linear in the vicinity of rj/xo- But if i]/xo is chosen too small, 

then it becomes very difficult to obtain reliable numerical estimates of the differences in T c . 

We found that rj/xo = 4 x 10~ 4 t was a good compromise in the nearly antiferromagnetic 

case. A smaller value of rj/xo was necessary for some of the nearly ferromagnetic results. In 

particular, we carried out the nearly ferromagnetic calculations for k 2 = 0.5 and g 2 Xo/t — 

10,5 and for g 2 Xo/t = 30 and k 2 = 2,3,4 with i]/xo = 5 x 10~ 5 t. We estimate that the 

accuracy with which we calculated the derivative dT c /d?7 in all cases is of the order of a few 

per cent. 

As in Reference 8 , we define a crossover frequency oo cross (q*) by AT c (q*, uJ cross ) = 0. 
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Similarly, one can define an optimal frequency cj opt (q*) as the frequency where AT C is 
maximum. uj opt (q*) then indicates which paramagnon frequency and wavevector q* con- 
tribute most to pairing. 

In order to make a comparison with the corresponding electron-phonon problem it is 
instructive to define a mass renormalization parameter \ z and interaction parameter A a- 
We define 

Xz = gy<^^(p~P / ^) >Fg(p,po (2Q) 

< 1 >FS |>: 

A jg v < M(p - p', "MpMpO >f5(p,pq (21) 

where 

V z (q,u) = g 2 x(q,cu) (22) 

and 

V p (q,u) = -^xM (23) 

?7(p) = sin(p x ) (24) 

for p-wave spin-triplet pairing (A = p) while 

V d (q,uj) = g 2 x(q,uj) (25) 
r}(p) = cos(p x ) - cos(p y ) (26) 

in the case of d-wave spin-singlet pairing (A = d). In carrying out the frequency integrations, 
we omitted the cutoff which was used in Eq. (10). The approximate effect of the cutoff is 
to multiply Aa and X z by a common factor which is weakly dependent on wavevector and 
k 2 . This can be ignored since we will only be interested in the quotient Xa/Xz- The Fermi- 
surface averages are given by 

<-'->F5(p) =/ 4^--- 5 ( e p-^ ( 27 ) 



(27T) 

d 2 p d 2 p' 
(2tt) 2 (2tt) 



< • • • >FS(p,p') = / , ^ 2 x 2 ■■■ 6 ( e P ~/Wp' -N ( 28 ) 
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In practice, we compute the Fermi-surface average with a discrete set of wavevectors on a 
lattice and we replace the delta function by a finite-temperature expression 

A&^S (29) 

5(e p - im) — > I/ P (l - / p ) (30) 

where / p is the Fermi function and iV is the number of wavevectors in The Brillouin Zone. 
Note that ±/ p (l - / p ) -> 5(e p - fi) as T -> 0. We used T = 0.1* and N = 128 2 . The 
finite temperature effectively means that van Hove singularities will be smeared out. The 
Fermi-surface average that appears in X z , Eq. (20), plays a role similar to that of a 2 F(uo)/uo 
in the case of phonon-mediated superconductivity. 



III. RESULTS 

The model consists of the parameters g 2 Xo/t, T s f/t, k and k. It is found experimentally 
that T s fK,Q is approximately constant. We shall use this relation to eliminate one parameter 
from the set and pick a representative value of the product T s /Kq. We put T s f = |t which 
corresponds to a temperature of 1000 K if the bandwidth is 1 eV, and n^a 2 = 12 for 
comparison with earlier work carried out on this model 2 . To obtain a representative value 
for the dimensionless coupling parameter g 2 Xo/t, we make use of the Stoner criterion. In the 
vicinity of the magnetic instability, <?Xo(Q, 0) ~ 1, where Xo(q, is the usual tight-binding 
Lindhard susceptibility and Q the ordering vector. With Xo(Q, 0) ~ N(0) ~ ^ where N(0) 
is the single-particle density of states at the Fermi level, the value of the coupling parameter 
g 2 Xo/t obtained is about 10. We stress that this is only an order of magnitude estimate of 
the coupling parameter. 

Figures 1 and 2 show AT c (q*,wo), Eq.(18), versus frequency uj for several values of (a) 
g 2 Xo/t, (b) k 2 and (c) q*. Figure 1 shows the results in the case of a nearly antiferromagnetic 
metal and Figure 2 shows the results in the case of a nearly ferromagnetic metal. By 
contrast to the phonon case where AT C is always positive, note that there are wavevectors 
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q* for which AT C is always negative, which means that spectral weight at these wavevectors 
is deleterious to superconductivity. And if spectral weight is added at sufficiently small 
frequencies uq, AT c is negative for all wavevectors q*. Finally, note that the overall scale 
of AT C for nearly ferromagnetic systems is larger than that for nearly antiferromagnetic 
systems. This indicates a larger sensitivity of the relative changes in T c with respect to 
changes in the spin-fluctuation spectral weight for nearly ferromagnetic systems. 
Figure 3 is a representative plot of 



as a function of the wavevector q* at which spectral weight is added for several values of u /t 
and for a nearly antiferromagnetic metal. AT C + indicates the region of momentum and fre- 
quency space where adding a small amount of spectral weight enhances the superconducting 
transition temperature. Figure 4 is a complementary plot for 



which indicates the region of momentum and frequency space where adding a small amount 
of spectral weight suppresses the superconducting transition temperature. These results 
should be compared with the corresponding plots, shown in Figures 5 (AT C + ) and 6 (AT~), 
for the nearly ferromagnetic case. We have chosen the values of u for graphs (a-c) to be less 
than uj opt (which, for these parameter values, is about 2.5£ in the nearly antiferromagnetic 
case and about 0.32t in the nearly ferromagnetic case) and the values of u in graph (d) to 
be greater than u opt . The region in wavevector space where AT C > grows monotonically 
with cu , starting from zero when u)q < Across, the crossover frequency corresponding to 
the incipient ordering wavevector. However the maximum value of AT C passes through 
an extremum at u = u opt , the optimal frequency corresponding to the incipient ordering 
wavevector. It is interesting to find that the size of the region AT C > is approximately the 
same in the nearly antiferromagnetic and nearly ferromagnetic cases if ujq is chosen close to 
the corresponding optimal frequency u opt , and that it is a significant fraction of the total 




(31) 




(32) 
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area of the Brillouin Zone. This suggests that except for very small values of k 2 , the pairing 
occurs mostly through the short-ranged magnetic fluctuations. Comparing the values of 
AT C in the frequency and wavevector scans suggests that removing spectral weight from low 
frequencies may be an effective way to enhance T c . 

We show in Figure 7 plots of AT C + for various values of k 2 in the nearly ferromagnetic 
case. The surfaces are very similar to one another in the small k 2 limit. This can also be 
seen in the AT C versus ojq curves. However, since our mean-field calculations are likely not 
accurate very close to the critical point for magnetic order, we cannot state categorically 
what the limiting values of AT C are. The AT C = contour in wavevector space is quite 
insensitive to the value of k 2 , but is, however, very sensitive to the value of u . 

Figure 8 shows how the crossover frequency u cross (q*) depends on the coupling constant 
g 2 Xo/t an d the correlation wavevector k 2 , for the optimum wavevector, namely q* = (tt,7t) 
for nearly antiferro magnetic systems and q* = (0, 0) for nearly ferromagnetic systems. The 
figure shows that oj ctoss (q*) when scaled by T c , is more robust to changes in g 2 Xo/t an d k 2 
in the nearly antiferromagnetic case than in the nearly ferromagnetic case. The value of 
u cross in the nearly antiferromagnetic case is about 3T C which is quite close to the value 2T C 
obtained in the Hubbard model 6 . The value of u cross in the nearly ferromagnetic case lies 
between 5T C and 9T C . It is interesting to observe that in the limit of small coupling and far 
away from the magnetic instability, oj cross /T c seems to approach a value common to both 
the nearly antiferromagnetic and nearly ferromagnetic cases. Figure 8(a) indicates that, for 
k 2 = 0.50 in the nearly ferromagnetic case, u cross /T c seems to depend on T s f and only 
through the product T s //tg, which is found experimentally to be approximately constant. 
Therefore u cross is roughly proportional to T c . But, for the same parameters, T c scales 
approximately linearly with T s 2 . Hence, for k 2 = 0.50 in the nearly ferromagnetic case, 
Across scales approximately linearly with T s f. 

One would like to understand the nature of the crossover frequency u cross in terms of 
phonon-problem-like parameters. Millis and coworkers 8 showed that, in the antiferromag- 
netic magnetic fluctuation case, 
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Across /T c ~ e^, (33) 

where 7^ is a d-wave effective interaction constant Id = j^- They assumed that the dy- 
namical magnetic susceptibility can be seperated into a wavevector-dependent factor and a 
frequency-dependent factor and they adopted an Einstein model to describe the frequency- 
dependent factor. They further assumed that the Eliashberg equations can be reduced to 
equations involving Fermi-surface quantities only. A similar calculation can be carried out 
in the nearly ferromagnetic case yielding 

Across /T c ~ e^, (34) 

with 7 P = yj. The results of the more complete calculation described in Section II are 

compared with these expressions for oJ cross /T c in Figure 9. Since 7a (A = p, d) is not an 

independent parameter in our model, we calculated log(a; cross /T c ) and I/7A for various values 

of k 2 . One sees that the expressions given in Eqs. (33) and (34) are a poor approximation to 

the value of oj cross /T c calculated from the Eliashberg equations. The dependence of u cross /T c 

1 

on 7a is much weaker than e 7 ^ and is not even monotonic in the nearly antiferromagnetic 
case. This finding is consistent with the inability to obtain a simple analytic expression 
similar to that proposed by McMillan to represent the T c calculated numerically via the 
Eliashberg equations 2 ' 4 . 

We show in Figure 10 the variation of the normalised optimal frequency u opt (q*)/T c with 
g 2 Xo/t an d k 2 in the nearly antiferromagnetic case, with q* = (tt, it), and nearly ferromag- 
netic case with q* = (0,0). In the nearly antiferromagnetic case u opt is approximately a 
constant equal to 10T C , which is close to the value it takes in the Hubbard model 6 . It is also 
of the same order of magnitude as 27rT c , the value it takes in the electron-phonon problem 5 . 
In the nearly ferromagnetic case u opt lies between 13T C and 22T C . 

As shown by Bergmann and Rainer for the phonon problem 5 , it is instructive to compare 
the frequency dependence of the spectral function with that of AT C . The results for the 
nearly antiferromagnetic case and q* = (tt, it) are shown in Figure 11 for various values of 

K 2 . 
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Figure 11 shows that in the small k 2 limit (strong coupling) the peak in the spectral 
weight lies below the optimum frequency oo opt , while for large k 2 (weak coupling) it lies 
above u) opt . This is analogous to the results of Bergmann and Rainer who found that the 
transverse phonon mode was below the optimum frequency u opt = 2ixT c for strong coupling 
superconductors such as Hg and above u opt for weak coupling superconductors such as In. 
The results of Figure 11 suggest an explanation for a maximum T c as a function of k 2 . As 
k 2 — > 0, while the d-wave component of the pairing interaction increases, the frequency 
at which the spin-fluctuation spectral weight is maximum becomes smaller than u cross , the 
frequency below which addition of spectral weight produces a suppression of T c . 

One might have expected that the n 2 for which T c is maximum would be such that the 
peak in the spin- fluctuation spectral weight coincides with uj opt . But for g 2 Xo/t = 10, the 
maximum T c is found at k 2 ~ 0.35, while the match between the peak in the spectral weight 
and uj opt occurs for k 2 ~ 1. The shift in the peak of ImV^ relative to u opt allows one to 
understand the trends in T c , but the arguments remain qualitative. 

Since the location of the peak in the spectral function does not depend on the value 
of g 2 Xo/t and our results show that the frequency u opt is approximately constant in the 
nearly antiferromagnetic case, one does not gain new insights by looking at the frequency 
dependences of ImVd(q*, u;o) and AT c (q*, aj ) for different coupling constants. 

In the nearly ferromagnetic case, since Im%(q — > 0, uo) = for uo ^ 0, one sees that the 
situation is quite different and thus the shape of the curves ImV^(q*, uo ) and AT c (q*,w ) 
for the optimum wavevector q* = will not be similar, unlike the phonon and nearly 
antiferromagnetic cases. 

IV. DISCUSSION 

Our results demonstrate a number of significant differences as well as similarities between 
the magnetic fluctuation- and phonon-pairing mechanisms. In the latter case, the interaction 
is local in space but non-local in time. This means that there should be almost no variation 
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of AT c (q*,c<Jo) with q* and therefore one only need worry about the dependence of AT C 
on frequency. In the case of magnetic pairing however, the interaction is non-local in both 
space and time and as a result AT c (q*,wo) exhibits a dependence on both wavevector q* 
and frequency c^o- 

In the phonon case 5 , addition of an infinitesimal amount of spectral weight at any non- 
zero frequency c^o results in an enhancement of T c , with an optimum frequency of 2ttT c . 
In the limit ooq — > 0, AT C — > 0, which is consistent with Anderson's theorem 9 for isotropic 
superconductors. On the other hand, for the magnetic interaction model, addition of spec- 
tral weight at sufficiently low frequencies results in negative values of AT C . The intuitive 
understanding of this result 8 is also related to Anderson's theorem. Very low frequency 
paramagnons effectively act as a static non-magnetic impurity potential that scatters the 
quasiparticles. For anisotropic superconductors, Anderson's theorem does not apply and the 
presence of non-magnetic impurities leads to a supression of T c . Also note that for certain 
wavevectors q*, AT c (q*,cj ) is always negative, regardless of u . 

There are also significant differences between the nearly antiferromagnetic and nearly 
ferromagnetic cases. For instance, the crossover frequency u) cross at which AT c (q*, aj ) = 
depends more weakly on the model parameters k 2 and g 2 Xo/t m the nearly antiferromagnetic 
case than in the nearly ferromagnetic case. (u cross is strongly dependent on q* in both 
cases.) Similarly, the optimum frequency u opt at which AT c (q*,w ) is maximum is more 
weakly dependent on the model parameters k 2 and g 2 Xo/t m the nearly antiferromagnetic 
metals than in the nearly ferromagnetic metals. Moreover, the overall scale of AT C is much 
larger in the former case than in the latter, indicating a much greater relative sensitivity of 
T c to changes in the spectral weight for nearly ferromagnetic systems. 

Our results also show similarities between magnetically mediated and phonon-mediated 
superconductivity. In all cases, the shape of the AT C curves are monotonic with a single opti- 
mum frequency (for each wavevector) and, in the high-frequency regime and for wavevectors 
such that AT C > 0, the curves behave similarly . The decrease of AT C with increasing ujq 
approximately mirrors the change in the interaction parameters A a as V e ff — > V e ff + SV e ff. 
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The relative positions of the peak u max in the spectral weight ImVd(q* = (7r, tt),u) and 
optimum frequency u opt in the strong coupling limit where 0J max < LV opt and weak coupling 
limit in which 0J max > ^o P t is similar in some respects to the findings of Bergmann and 
Rainer 5 for the phonon problem. In that case, the transverse phonon frequency is below the 
optimum one for strong coupling superconductors and above in the weak coupling case. 

Figure 2(c) shows a feature that may partially explain the difficulty in observing super- 
conductivity on the border of itinerant ferromagnetism. Our results (see Figure 2(c)) show 
that adding spectral weight at large wavevectors is detrimental to T c . In other words, an- 
tiferromagnetic correlations act to suppress pairing in nearly ferromagnetic systems. In the 
presence of a lattice, one would generically expect some enhancement of the magnetic re- 
sponse at large wavevectors due to nesting features in the Fermi surface for nearly-half-filled 
electronic bands. As explained in References 2 ' 4 , magnetic pairing on the border of long- 
range ferromagnetic order is not as robust as magnetic pairing in nearly antiferromagnetic 
systems. The antiferromagnetic correlations inherent to the presence of a crystal lattice 
tend to suppress magnetic pairing in nearly ferromagnetic systems, thus making it even 
more difficult to observe. 

V. OUTLOOK 

The dependence of the superconducting transition temperature on model parameters and 
the role played by dimensionality for Ornstein-Zernike like spin-fluctuation spectra has been 
studied in some detail 2,4 . However, these calculations reflect an integral property of the 
dynamical susceptibility x(q, a;). In this paper, we presented a study of how different regions 
in wavevector q and frequency cu individually contribute to T c , giving us novel insights into 
the magnetic interaction model. 

The calculations were carried out by adding infinitesimal amounts of spectral weight 
at different wavevectors and frequency and determining the resulting changes in T c from 
the numerical solution of the Eliashberg equations. We stress that the work reported here 
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provides insight only into how different regions in q and uj contribute to pairing for the 
particular magnetic-fluctuation spectrum x(q, cu), Eq. (10). One may infer certain trends 
in T c for small changes in the spin-fluctuation spectrum, but it may not be warranted to 
extrapolate our results to large changes in the magnetic spectrum. Also, the calculations 
are carried out at T c . They may give one some idea of how the robustness of pairing is 
affected as one goes below T c and the electronic gap induces changes in the momentum and 
frequency structure of x(q, but this involes another extrapolation to large changes in the 
magnetic spectrum which may not be warranted. 

We found that addition of spectral weight at or near the incipient ordering wavevector 
in the nearly antiferromagnetic and nearly ferromagnetic cases lead to an enhancement 
of T c , except at low frequencies where it leads to a suppression of the superconducting 
transition temperature. However, addition of spectral weight far away from the incipient 
ordering wavevector results in a lowering of T c regardless of frequency. These results are in 
stark contrast to those obtained for phonon-mediated superconductivity where addition of 
spectral weight at any non-zero frequency leads to an enhancement of the superconducting 
critical temperature. 

The theoretical framework presented here to describe systems on the border of magnetism 
can be translated to describe systems on the border of other types of instabilities, such as 
charge- density- wave or ferroelectric instabilities. The same type of calculations could be 
carried out in those other cases, and compared to the results of this paper, shedding light 
on the similarities and differences between the many possible pairing mechanisms. 
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FIGURES 

FIG. 1. Frequency dependence of AT C = f^^^ for various values of (a) the coupling parameter 
g 2 Xo/ti (b) the inverse correlation length re 2 and (c) the momentum transfer q* in the case of a 
nearly antiferromagnetic metal with a d x i_ y i superconducting state symmetry. The characteristic 
spin-fluctuation temperature is T af = 0.67 and reg = 12. 

FIG. 2. Frequency dependence of AT C = for various values of (a) the coupling parameter 

g 2 Xo/t, (b) the inverse correlation length k 2 and (c) the momentum transfer q* in the case of a 
nearly ferromagnetic metal with a p-wave superconducting state symmetry. The characteristic 
spin-fluctuation temperature is T s f = 0.67 and Kq = 12. 

FIG. 3. Wavenumber dependence of AT C + = min{^-^^,0} for various values of uJo/t for a 
nearly antiferromagnetic metal with a d x 2_ y 2 superconducting state symmetry. The other pertinent 
parameter values are g 2 Xo/t = 30 and k 2 = 0.5. 

FIG. 4. Wavenumber dependence of AT" = mm{^-^^,0} for various values of uo/t for a 
nearly antiferromagnetic metal with a d x 2_ y 2 superconducting state symmetry. The other pertinent 
parameter values are g 2 Xo/t = 30 and k 2 = 0.5. 

FIG. 5. Wavenumber dependence of AT+ = max{^^^,0} for various values of uo/t for a 
nearly ferromagnetic metal with a p-wave superconducting state symmetry. The other pertinent 
parameter values are g Xo/t = 30 and k 2 = 0.5. 

FIG. 6. Wavenumber dependence of AT" = minj^^^jO} for various values of uj^/t for a 
nearly ferromagnetic metal with a p-wave superconducting state symmetry. The other pertinent 
parameter values are g 2 xo/t = 30 and k 2 = 0.5. 

FIG. 7. Wavenumber dependence of AT+ = max{^^, 0} for various values of the dimension- 
less parameter n 2 for a nearly ferromagnetic metal with a p-wave superconducting state symmetry. 
The other pertinent parameter values are g 2 Xo/t = 30 and ujo/t = 2. 
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FIG. 8. The dependence of the crossover frequency uj cross , normalised by the superconducting 
transition temperature T c , on (a) the coupling parameter g 2 Xo/t and on (b) the inverse correlation 
length k 2 . Both graphs exhibit results obtained in the nearly antiferromagnetic (NA) case and 
for three values of the parameter pair (T 8 f, Kq), keeping the product T s fK,Q constant, in the nearly 
ferromagnetic (NF) case. Throughout we have put q* = (it, it) in the nearly antiferromagnetic case 
and q* = (0, 0) in the nearly ferromagnetic case. 

FIG. 9. Comparison of our results for the crossover frequency in (a) the nearly antiferromag- 

i 

netic case (circles) and (b) the nearly ferromagnetic case (circles) with the result cu cross /T c ~ e 1 ^ , 
Eq. (33) (squares). We have put the prefactor in this formula equal to one for the sake of com- 
parison. We have calculated oj ctoss with g 2 Xo/t = 10, although the effective interaction constant 
7A is independent of g 2 Xo/t- The data is parametrised left to right by increasing k 2 from (a) 0.10 
to 4.00 and (b) 0.10 to 1.00. Throughout we have put q* = {it, it) in the nearly antiferromagnetic 
case and q* = (0, 0) in the nearly ferromagnetic case. 

FIG. 10. The dependence of the optimal frequency uj opt , normalised by the superconducting 
transition temperature T c , on (a) g 2 Xo/t and (b) k 2 . We have put T s f = 0.67, Kq = 12 and 
q* = (-7T, it) in the nearly antiferromagnetic case and q* = (0, 0) in the nearly ferromagnetic case. 

FIG. 11. Comparison of AT C = jr^j- (solid line) and ImV^ (dashed line) for various values of 
the parameter k 2 for a nearly antiferromagnetic metal. The functions are plotted for a value of the 
coupling g 2 Xo/t = 10 and a wavevector q* = (it, it). 
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